Relaxation time ansatz and shear and bulk viscosities of gluon matter 



A. S. Khvorostukhin, 1,2 V. D. Toncev, 1 and D. N. Voskrescnsky 3, 4 

1 Joint Institute for Nuclear Research, 141980 Dubna, Russia 
2 Institute of Applied Physics, Moldova Academy of Science, MD-2028 Kishineu, Moldova 
s National Research Nuclear University "MEPhI", Kashirskoe sh. 31, Moscow 115409, Russia 
''GSI, Helmholtzzentrum fur Schwerionenforschung GmbH, Planckstrasse 1, 64291 Darmstadt, Germany 

Shear and bulk viscosity-to-entropy density ratios are calculated for pure gluon matter in a 
nonequilibrium mean-field quasiparticle approach within the relaxation time approximation. We 
study how different approximations used in the literature affect the results for the shear and bulk 
viscosities. Though the results for the shear viscosity turned out to be quite robust, all evaluations 
of the shear and bulk viscosities obtained in the framework of the relaxation time approximation 
can be considered only as rough estimations. 
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INTRODUCTION 



High-energy heavy-ion collisions at RHIC and LHC en- 
ergies have shown evidence for a new state of matter 
characterized by very low shear viscosity to entropy den- 
sity ratio, r\j s, similar to a nearly ideal fluid (see [l|, Q)- 
Lattice calculations indicate that the crossover region be- 
tween hadron and quark-gluon matter has been reached 
in these experiments. On the other hand, for the pure 
gluon SU(3) theory lattice calculations demonstrate the 
occurrence of the first-order phase transition (see 0,11]). 
Recently, there appeared gluon lattice data on ratios of 
the shear 0, 0] and bulk 0, viscosity to the entropy 
density. 

Among various existing phcnomcnological approaches, 
quasiparticle models are used to reproduce results ob- 
tained in the lattice QCD (see 0-fHJ)- In the case of glu- 
odynamics, quasiparticle models rely on the assumption 
that for a temperature T below the critical one, T < T c , 
the system consists of a gas of massive glueballs and, for 
T > T c , the system consists of deconfincd gluons. Pcrtur- 
bative estimates of the shear and bulk viscosities in gluo- 

and nonperturba- 
within the quasi- 
particle models with temperature-dependent masses in 
the relaxation time approximation with some additional 
ansatze resulting in essentially different expressions for 
the bulk viscosities. 



dynamics were performed in 0, 13, 
tive evaluations were made in 15 



II. THE QUASIPARTICLE MODEL 

We start with the expression for the gluon energy- 
momentum density tensor in a relativistic mean-field 
model, 



= / rfr 



f(t,r;E,p)+g^B[m 2 (p)} , (1) 



obeying the conservation law 



d^ v = 0. 



(2) 



Here 



dY 



d 3 p 

'(27T) 5 



, p^ = {E[m 2 (ip)],p} 



for gluons with a degeneracy factor d g = 2(N 2 — 1) = 16, 
N c = 3; ip is a mean field; summation over the repeated 
indices is implied. In this model the quasiparticle energy 
is given by 



E = \/p 2 + m?{p) 



(3) 



We assume that the quasiparticle distribution function / 
obeys the kinetic equation 

E-tydpf-VEVtf = St/, (4) 

where St/ is the collision term satisfying the condition 

dTE St/ = 0. (5) 



Here we continue to exploit the mean-field-based quasi- 
particle model with parameters fitted in [l6| to fulfill the 
modern lattice data. Using two possible ansatze for the 
collision term in the relaxation time approximation we 
derive expressions for the shear and bulk viscosities to 
entropy density ratios for the pure gluon SU(3) theory 
at T > T c and compare the results with the lattice data. 
Our results demonstrate ambiguities in calculation of the 
bulk viscosity within the relaxation time approach. 



Applying ([2} for v = and using (fTJ) and §5§ we derive 
the consistency condition for the effective bag constant 



dB 

dtp 



where 



dm 
dip ' 



E J 



(6) 



(7) 
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is the scalar density. Note that in these expressions the 
quasiparticle energy is a functional of the distribution 
function E = E[f]. In thermal equilibrium relation ([6]) 
coincides with the thermodynamical consistency condi- 
tion derived in and used in the model (l6j | . 

The local equilibrium distribution function for a gluon 

is 



l.eq. 



(8) 



Here rf eq . = (£ Leq - 



P), 



E le< i- = Jp 2 



m 2 {<p Le( i[T(t,r)}}, 



and the four-velocity of the frame is ~ [l,u(t, r)] for 
|u| -C 1. Following [la], we use the pole mass 



(m Leq -) 2 =mV' eq '( T )] 
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„2 



fl 2 (T) T 



(9) 



as the gluon quasiparticle mass in local equilibrium. 
The temperature-dependent strong interaction coupling 
in the next to leading order is given by 



9 2 (T) 



48^ 2 



.(10) 



All previous papers exploited g 2 up to the leading or- 
der using only the first term in the denominator of (|10p . 
In Ref. [16l |. working within this leading order we fitted 
the parameters A and T s to fulfill the new lattice data 
[H for the reduced pressure, energy, enthalpy and trace 
anomaly. However, we found that in the perturbativc 
limit [at In ln( x< '- T T Ts ' > ) 2 ^> 1] the correction including 
the double-logarithmic term yields a larger contribution 
to dm 2 /dT 2 than the leading logarithmic term. Since 
one of our aims in the given paper is to improve our fit 
in the perturbativc limit, we keep the double- logarithmic 
correction and tune parameters to reproduce thermody- 
namic characteristics of the system. Our fit [l6| with g 2 
calculated up to leading order was performed with the 
parameters T s /T c = 0.5853, T c = 265 MeV, and A = 3.3. 
The new fit with g 2 calculated up to sublcading order 
yields T s /T c = 0.015, T c = 265 MeV, and A = 1.53. 
Our previous fit of the reduced pressure, energy, enthalpy 
and trace anomaly (i~6| proves to be almost unchanged. 
Therefore, we do not redraw those figures here. 

Generally, the running coupling constant a(T,r) = 
g 2 (T,r)/(4ir) depends on the temperature and distance 
and differs from the effective one determined as a(T) = 
5 2 (T)/(47r) following Eq. ([TO]). The quadratic rise of 
a(T = 0,r) (see solid line in Fig. 6 of [201 ]) is a nonper- 
turbativc effect that stems from the linear rising string 
tension term in the potential, while at small distances 
the logarithmic weakening of the coupling is visible and 
at sufficiently small distances it reaches the perturba- 
tive behavior (asymptotic freedom). At finite tempera- 
tures, a(T, r) follows this zero-temperature behavior to 




FIG. 1: (Color online) The effective running coupling a(T) = 
g 2 (T)/ (4-7r) as a function of T/T c . The points are the lattice 
data [20] (see explanation in the text). The solid line is the 
result of our new fit (|10|l . The dash-dotted line is our (leading- 
order) fit [l|J. The dashed curve is the (leading-order) fit of 
Ref. [!l with T s /T c = 0.73, T c = 271 MeV, and A = 4.3. 



relatively large distances, before Debye screening sets in, 
leading to a maximum and a decrease at larger distances 
(see Fig. 6 of Ref. [I(J). The value a is a temperature- 
dependent coupling constant determined with the help 
of the Debye screened potential at large distances [2fJ] . 
Thus, the behavior at large distances of the coupling, 
which one may use to calculate thermodynamical charac- 
teristics in quasiparticle models, is not defined uniquely. 
Therefore, in Fig. [1] along with the effective running cou- 
pling constant a (circles) the coupling constant a max 
(squares) is plotted to be defined by the maximum in 
a(T,r) [2fJ. As can be seen, both estimates of the run- 
ning constant noticeably differ near the critical temper- 
ature and tend to coincide at high temperatures. In Fig. 
[TJ we demonstrate also that the running constant a(T) 
calculated up to the next to leading order with the help 
of Eq. (ITU1) in the given work (solid line) and those cal- 
culated up to the leading order with the parameters from 
flr| (dash-dotted line) and from [l?} (dashed line) devi- 
ate only little from each other, describing reasonably the 
lattice data of Ref. [2(| ■ Note that the parametrization 
from Ref. 



17| is fitted to the old lattice data on thermo- 



dynamical quantities and does not describe properly the 
new lattice data B. 



III. SHEAR AND BULK VISCOSITIES 

We define the shear, 77, and bulk, £, viscosities as 
coefficients, entering into the variation of the energy- 
momentum tensor in the local rest frame: 



ATij = -C SijV -u + rjWij, 

2 _ 

with Wij = diUj + djUi + -Sij V ■ u 

o 



(11) 
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Here Latin indices 1, 2, 3 correspond to the spatial com- 
ponents. Summation over the repeated spatial indices is 
implied. To find the shear viscosity, 77, we set i += j in 
(fTTj) . To find the bulk viscosity, £, we substitute i = j in 
([IT]) and use the fact that I^ oq = 3P lcq . 

The operation A in (fTTj) needs explanation. According 



21 



in quasi- 



to the work of Abrikosov and Khalatnikov 
particle Fermi liquid theory one usually exploits the fact 
that the following combination enters into the original 
Landau collision term: 

6 (^> a[/] ) ' + /(W + f&i)] 

-f(E 3 )f(E 4 )[l + f(E 1 )][l + f(E 2 )}} , (12) 

where E a are the functionals of the exact nonequilibrium 
distribution function f a and a = 1,2,3,4. The term in 
curly braces term is zero for f l ' Gq (E a ) given by 



/'■ cq '(£) 



with?/ 1 = [E(p),p\ , (13) 



where E(p) is determined by Eqs. ([3]) and ([6]) and 
depends on the nonequilibrium distribution. Thus, 
Stf l - cq -(E) = 0. However, following [lf| the same col- 
lision term should vanish also for local equilibrium, i.e., 
St/ 1 ' eq (£ ,Loq ) = 0. Let us demonstrate this by an exam- 
ple of the 2—^2 processes, as they are described by the 
four-momentum kinetic Kadanoff-Baym equation. Then 
the local collision term renders u3 



c ioc. = d 2 



d A Pi d 4 p 2 d 4 p 3 2 



(2tt) 4 (2tt) 4 (2tt) 4 1 
xS 4 (p + pi -p 2 - p 3 ) A(x,p)A(x,p 1 )A(x,p 2 )A(x,p 3 ) 
x {f(x,p 2 )f(x,p 3 )[l + f(x, Pl )][l + f{x,p)} 
-[1 + f(x, P2 )][l + f{x,p 3 )]f{x, Pl )f(x,p)} , (14) 

where M is the properly normalized matrix element and 
A = —2ImG R is the spectral function (density of states) 
that implicitly depends on the corresponding distribution 
function /, G R is the retarded Green function, andpo and 
p are independent variables. The Landau quasiparticle 
collision term is obtained from here provided one sets 
A[f] = 2nS[p 2 -p 2 -m 2 - RcX R (p, [/])], where m is 
the free particle mass and is the retarded self-energy. 
Here the local equilibrium distribution fulfills the relation 



/ LCq -(.T, P + (Z )[l + / 1 - Cq -(x,p)] 

= [f Lcq -(x,p) - f- cq -(x,p + q)}f- cq -(x,q), 



where 



f cq (x,p) 



exp 



T(t,r) 



(15) 



(16) 



i.e., in contrast to © and ([TB"]) . po and p are indepen- 
dent variables. With the help of this relation we can see 



that the term in the curly braces in (|14[) is zero indepen- 
dently of values of A. Thus, C loc ^ = for two distribu- 
tions: F = iG- + = A[/ 1 - eq ']/ Le<1 ' and F = A[f]f heq ' 
(whereas the real solution of the kinetic equation is 
F = A[f]f). In our model m does not depend on p, and 
the Landau quasiparticle collision term is obtained pro- 
vided one sets in A[f- cq -} = (TT/E Lcq -)S(p - E lcq ) 
(or A[f] = / E[f])8{p Q - E[f})). Note that in the Lan- 
dau quasiparticle kinetics one neglects memory effects; 
therefore, we indeed may use the local approximation in 
Eq. d). 

Then, returning to our quasiparticle model we can ex- 
pand the distribution function near f^^i-^E 1 ^-), i.e. (8]), 
performing variations as 



Sf = f(E) -/ Lcq -(£ Lcq -), 

ST^ = T^ v [f(E)\ — T^ v [f heq (E Leq ')] 



(17) 



i.e., with A = (5 in Eq. (jlip . and alternatively we can 
introduce the operation A as A = 8, where 

Sf = f(E[f])-f l - cq -(E[f]), (18) 
6T" V = T^[f{E)]-T^[f Lcq (E)}, 

where / hcq ' (£[/]) is given by Eq. (13]). Subtracting (17]) 
from (fT8l) we obtain the relation 



Of 



= ^) SE. 



l.eq. 



(19) 



Below we use both possible choices of the A operation. 
Performing variations in (TJ we find 



(20) 



where we used the fact that after integration over angles 
/ dTpiPj(f /E 2 )AE = 0; the diagonal terms 

AT U = J dT^Af - J dT^fAE - SAB, (21) 

where in the case A = 8 the_second and the third terms 
vanish since B = B(E) and SE = 0; and 



AT 



00 



dTEAf. 



(22) 



To derive (|2"2"]) , we used the consistency relation (5]). 

Since in the local rest frame the viscosities are pre- 
factors at small variations of the velocity, in the linear 
approximation used their values should not depend on 
how variations were performed (using operation 8 or 8), 
provided all derivations are made without additional as- 
sumptions. These values are usually computed with the 
help of the local equilibrium distribution functions (|8]). 
However, in practice one uses additional approximations, 
e.g., great simplification arises if one exploits the so- 
called relaxation time approximation or more precisely 



4 



the relaxation time ansatz. Moreover, one can impose 
the so-called Landau-Lifshitz or somewhat different con- 
dition. Below we show that the use of different ansatzc 
can bring significant differences in the resulting values of 
the kinetic coefficients derived by means of the operations 
6 and 5. 

In our previous work (l6| . to find shear and bulk vis- 
cosities of the gluon and glueball matter we considered 
two relaxation time ansatze (see Ref. [23[) choosing the 
collision term as 



st/ = -<y//r, 

and in a different form 

St/ = -J//f, 



(23) 



(24) 



where r and r are, in general, different values and in gen- 
eral they are energy- and momentum-dependent quanti- 
ties. All previous works SEMI used averaged values 
for the gluon relaxation time. In the latter case, as we 
show below, the consideration can be performed consis- 
tently whereas for the energy- and momentum-dependent 
relaxation time, problems can arise [e.g., the two condi- 
tions (|4"5l) might be not fulfilled simultaneously]. Thus, 
below we set r and r to be constant unless told otherwise. 

Additionally, for relativistic systems one usually uses 
the approach of Landau and Lifshitz where u M is defined 
as the four- velocity of the energy transport. Thus we 
require 



AT M = 0. 



(25) 



Then in the local rest frame the energy should satisfy the 
Landau-Lifshitz condition 



ST 



oo 



0. 



(26) 



Also, performing all calculations at fixed exact particle 
energy, i.e., with the help of the (5-operation, we can re- 
quire fulfillment of the condition 



ST m = 0. 



(27) 



Thus, below we study three possibilities: 
ansatz I: when the right-hand side (r.h.s.) of the ki- 
netic equation is presented in the form (|23[) and the 
Landau-Lifshitz condition ([2l)|) is imposed, 

ansatz II: when the r.h.s. of the kinetic equation is pre- 
sented in the form (f2~4"|) and the condition (|27[) is imposed, 
and 

ansatz III: when the r.h.s. of the kinetic equation is 
presented in the form (f2"4")l but the Landau-Lifshitz con- 
dition (|26| is imposed. 



A. Relaxation time ansatz I 

Let us first consider the kinetic equation with the 
r.h.s. in the form (|23]l and use the Landau-Lifshitz con- 
dition 4H (ansatz I). Replacing / ~ /i-eq-(£i-eq.) [j. e-) 



at \Sf \ < /Leq-(£ji-eq.)] j on the left-hand side (l.h.s.) of 
the kinetic equation ((4]) wc find 



l.cq. 



(28) 



/(! + /) 



l.eq. 



PkPiWki 



, /(! + /) 
T 



2TE 

Introducing two terms 

Sf = Sf[W M ] + Sf[V-u] 
and using ((23)) for the r.h.s. of ([4]), we get 



l.eq. 



./(! + /) 



2TE 



l.eq. 



PkPlWkL, 



and for Sf[V • u] we find 



J/[V • u] 



, /(! + /) 
T 



l.eq. 



V • u. 



V • u. 



(29) 



(30) 



(31) 



Here the equation-of-state-dependent Q factor (see [Tl|, 

in Hi) 

is given by 



Q 



3E 



E-T 



dE 
df 



and 



3<9T 00 



(32) 



(33) 



is the speed of sound squared of the local equilibrium 
system. To get (|3Tj) . we used 



(34) 



9£ ;l.eq. ^ dE \.e q . ^l.cq. Q T 

dip lcc i- dT ~dt 



dt 



Substituting ([50)) into ([2"U)) and comparing with (TTTJ), we 
easily find the expression for the shear viscosity, 

where for convenience we introduced the notation 

(<f(p)) = j dT [$(p)/(l + /)] Lcq ' . (36) 
with respect to the mean field we 

l.eq 



Varying Eq. 
obtain 



5 ip 



A 



dm 

dtp 



l.eq. 



with 



A _ A = ( dm 

dip 



l.eq 



E \.cq. 



(37) 



,2 \ l.eq 

d T~ 



dip 2 



l.eq. 



f d 2 B 
\d^ 



l.eq. 



(38) 
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Calculating d 2 B/dip 2 with the help of (|6|) we get 



T \d(p 



m dm ( dm\" I m 2 



•(39) 



Making use of 5m = (dm/dif)5ip from Eqs. (|37p and ([55)1 
we have 



dm 



rp2dm f r ly-§± 
1 dT J E 



Then from Eqs. (J37J and (gTJ) we find 

' p 2 - K 



5T V , 



dT 



E l.cq. 



Sf, 



(40) 



(41) 



where 



K = -Tm 1 



dm 
~dT 



P 

E 2 



/i\ rj^dm I m 2 



(42) 



Using ([TT]) we obtain the final expression for the bulk 
viscosity: 



s T\3E yy 



Note that for r 



= const the variable shift 
Qt^Qt + bE, 



(43) 



(44) 



suggested in [23( and then used in [13, LL8( to satisfy the 
Landau-Lifshitz condition (|26[) . is not required here, since 
Eq. ((261) and the condition / dT E St f = 0, i.e., 



(tQE) = and (QE) = . 



(45) 



are fulfilled with Q given by Eq. (|32j) . This statement is 
easily checked with the help of the relation 



(c l.e qf 



(P 2 ) 



3 {E 2 } 



Tm— 



(1) 



(46) 



To reproduce the second relation (|45l) we used the result 
([3T|) and Eq. (f23]> . 

Moreover, in the case of a momentum-dependent relax- 
ation time and also in the case of many particle species 
with different values of the relaxation time, r a , the re- 
placement Q a T a — > QaTa + bE a does not generate new 
solutions St f a = 0, provided the collision term is pre- 
sented within the relaxation time approximation. Thus, 
we note that the replacement Q a r a —¥ Q a T a + bE a gen- 
erates new solutions St/ a = and thereby the shift is 
meaningful only if the exact form of the local collision 
term St/ is used. 

Nevertheless, compared to the results of other works, 
performing the shift (|4"4"]l one can present Eq. (|4"3"]l in the 
form 



< = -f(g(? 2 -*-^ 2 ) 



(47) 



with 
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(P 2 )-K(l) 
HE 2 ) 



(48) 



where following (|4"5j) the last term in Eq. ([4"T]) is zero pro- 
vided r is constant. We note that, although the variable 
shift (|44p has been performed, Eq. (|4Tt does not yield 
the quadratic form exploited in Refs. jlfl [l8l . [23| . Only 
for m = const do we recover the quadratic form for the 
bulk viscosity 24| . We also note that for a temperature- 
independent g 2 , expression (|47[) would yield C = 0. This 
circumstance gives additional justification to our inclu- 
sion of the double- logarithmic correction to g 2 (T) in Eq. 
(|T0)) in order to correctly take into account the m(T) 
dependence. 

Now let us prove that expression (|4T|) is nonnegative. 
For that we rewrite Eq. (|4T|) as 



T 

9T 

x 



3c 2 T 2 



dm 
dT 2 



K 



1 

E 2 



(l) 2 



From here for m = const we obtain 



C(m = const) 
From Eq. 



9T 



1 

E 2 



OH 

(E 2 ) 



(49) 



(50) 



that m 2 > T 2 % 



46| and the condition 3c 2 < 1, it follows 
and K > — m 2 . Then making use of 
C(m = const) > we find that ( > 0. 

Note that from relations ([15]), ([3T]l and (gDJ it follows 
that the collision term ((23]) can be expressed as 

(51) 



St/[V • - -«5/[V ■ u]/t = 5f[V ■ u]/t{E) 
with the energy-dependent relaxation time 
Tm^{r%)/{rEQ) 



r(E) 



1 



(1) 



Tm—L-K 

'"'dT \ W 



(52) 



Expression (f5T|) (see the second equality) has the same 
form as being postulated in ansatz II. However, note 
that to satisfy the Landau-Lifshitz condition in the form 
(|26p . we assumed that t is an averaged value and the 
relaxation time t(E) determined by the second equality 
([ST1) proved to be a certain (momentum-dependent) func- 
tion depending on the constant value r. Thus, starting 
with ([2U|) for r = const, we arrived at ([Ml) with r = t(E). 
Also in passing we demonstrated that the collision term 
(f2"4"| [as well as (|23[)] vanishes in the local equilibrium 
state [i.e., for / = f^i^E 1 ^-)}. 



B. Relaxation time ansatz II 

Let us now assume that the collision term is given by 
Eq. (|24[) and perform all variations at a fixed nonequi- 
librium energy, i.e., we apply the S operation to all quan- 
tities. Then SE = and 5B(E) = and the distribution 



6 



function / l oq (i?) is given by (fT3j) . Since we need to keep 
only terms linear in diU^ on the l.h.s. of the kinetic equa- 
tion, we may use there <j8j> instead of (p~3]l . Then using 
Eq. U2SJ| for the l.h.s. of $D and Eq. for the r.h.s. 
we find 



and 



Sf[W k i] 



; /(! + /) 

2TE 



l.eq. 



PkPiWki 



, /(! + /) 
T 



l.eq. 



V • it. 



(53) 



(54) 



Applying here relation (|19j) and averaging the result over 
angles we arrive at 



5T k = ST lk , for i y£ k, 



and also 



ST io = ST 



(55) 



(56) 



By using the consistency condition, relation ([T9")) after 
partial integration gives 



(57) 



ST U = I dTp 2 Sf/E 



Exploiting ansatz II we replace the Landau-Lifshitz con- 
dition ([25)1 by condition ([2"7]) . The last condition differs 
from (|26|) by the terms linear in 8m oc Vit. We observe 
that the condition ([27)) with Sf satisfying ([54"[) is fulfilled 
provided r is assumed to be a momentum-independent 
(averaged) quantity since, as in ansatz I, Eq. (|27[) is re- 
duced to the same second condition (|45[) proved above 
for Q from Eq. (pjl. 

Then we immediately recover expressions for the shear 
viscosity, 



(58) 



and the bulk viscosity, 

c = 



1 ( tQp 2 
T \ 3E 



(59) 



For m = const, Eq. (|59|) tr ansforms to ([50]) and repro- 
duces the result of Ref . [24| . 

If for completeness of consideration, the shift 

Qt^Qt+Ie (60) 

is performed, Eq. (f59| can be presented in the form 



1 ItQp' 



3E 



Kim/iEi\. (61) 

T (tE 2 ) \ 3 ; 



The last term in Eq. (|6Tj) vanishes if one uses the aver- 
aged value of t. Again, although the variable shift ([6"0")) 
has been done, Eq. ([611) does not yield the quadratic form 
exploited in Refs. jl7l li~8l. [23j. Also, one should stress 
once more that the variable shift is meaningful only if the 
exact collision term without memory terms is used. 

Note that the result g>3) follows from flUJ) for K = 0. 
Thereby, the final result (|5T)]) is nonncgative. 

We assumed that values r and r entering into the col- 
lision terms ([23"]) and ([M|) are averaged quantities (con- 
stants) evaluated by means of the cross section. In this 
case, we cannot distinguish between r and r and should 
set t = t = t with f being an averaged relaxation time. 
Thus, with such a simplified approach we are unable to 
distinguish which expression, (|43[) or (|59[) . is more prefer- 
able. Then for r\ both approaches (with A = 6 and 6 = 6) 
result in the same expression. 

Exploiting ansatz I in [T(| we further used an addi- 
tional ansatz, namely, in performing variations we did 
not vary quantities which depend on the distribution 
function only implicitly, such as E and B(E). Doing 
so we arrived at Eq. ([S"9"]) . Thus, our results obtained 
in (l6| within ansatz I are actually equivalent to those 
derived with ansatz II in the given work. As follows from 
the comparison of Eqs. (|4"3")l and ([5^1 the approxima- 
tion used in [l6T | within ansatz I is appropriate only for 
\K\ -c m 2 . The latter inequality holds for nonrelativis- 
tic quasiparticles but it becomes a poor approximation 
in the relativistic limit. As follows from Eq. © and 
Fig. [TJ in our case the nonrelativistic approximation is 
not applicable. Therefore, calculations with one and the 
same value f performed following Eqs. (|4"3")) and ([5^| may 
differ significantly. 

Combining (jT^j) and ([3T))) . we find 



Sep 
L- 1 



dm 



l.eq. 



L m— / dT— : — 6f . 

dip J J E he * 



1 



(62) 



A" 1 



1 / dm\ m 



T \ dip J E 2 



l.eq. 



J dT m dm , 
[d^T 2 



and 



8m 



T 2 dm 



dT 



E 



(63) 



(64) 



Using now relations ([H)|) , and ([51[) we are able to express 
the collision term ([24]) in the form 



St/[V 



Sf[V 



with 



-(E) = 



1 - 



-(E) 



dm /rQ\ 
± dT \ E i 



-EQ(\) 



(65) 



(66) 
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Thus, if one assumes that Eq. ([241) holds for r ~ 
const, then one recovers Eq. (|65p with the momentum- 
dependent value of t(E). In contrast, if one assumes that 
Eq. ([2"B"]) holds for r ~ const, then one recovers Eq. (|5Tj) 
with the momentum-dependent value of t{E). 



C. Ansatz III 

It was accepted in Ref. [23| that STu = 
/ dTp 2 5f/E lc i- but the Landau-Lifshitz condit ion was 
used in the form ([2T]) . Moreover, for simplification the 
authors [23[ assumed a specific form of the nonequilib- 
rium distribution function / = e^ E ^^ T ^ (in the Boltz- 
mann limit /< 1) introducing variations of the tempera- 
ture ST. The mentioned procedure results in a quadratic 
form for the bulk viscosity. 

Let us demonstrate how one can derive the same 
quadratic form for ( but slightly more consistently. If 
we ignored the Landau-Lifshitz condition, then using the 
first line of (|57[) and the kinetic equation within ansatz 
II ([24]) . we would recover expression ([5§|) . On the other 
hand, using the second line of (|57|) . we can expand the 
energy-momentum tensor density near the local equilib- 
rium state and, therefore, exploit Eqs. and (|20|) — (f22|) 
for A = <5, taking, nevertheless, the collision term in the 
kinetic equation in the form ([24"]) . 

We express ([26]) through Sf using ([62]) . Then 

df m dm^ 



E dT J 



One can easily see that expression ([57]) is not zero if ([54"]) , 
i.e., the Landau-Lifshitz condition ([2l)]). is not fulfilled. 
Performing the variable shift (|60| one can find the pa- 
rameter b to satisfy condition (|26[) . Thus, we find 



: <2 [E 



Tm Am 
~E~dT 



(68) 



However, we observe that after the variable shift the ki- 
netic equation within the initially assumed Eq. (|24l) is no 
longer fulfilled. Nevertheless, we can find the new kinetic 
equation 



St/ 



5f + b—f L ^(l + f-^)V-u 



,(69) 



being in agreement with the performed variable shift and 
the Landau-Lifshitz condition. 

Obviously, the modification ([M]) compared to the re- 
sult ([24]) does not affect the result ([58]) for the shear 
viscosity. Combining ([57]) . ([6"9"]h and ([S5]) one arrives at 
the result of Ref. 23( for the bulk viscosity: 



< 



(rQ 2 ) 



(70) 



For m = const Eq. (|70l) transforms to ([50]) , reproducing 
the result of Ref. |24{ . 

Note that by means of Eq. ([28]) we can rewrite ([69]) as 



St/[V 



u = — 



with 



T 3 (E)=T 1 




(71) 



(72) 



Thus, the collision term ([M]) reaches zero for / = 
/i cq.(£;[V • u]) given by ([15]). 

Finally, let us estimate the difference between values (j 
[Eq. O] and Cm [Eq. ([70])] . After simple but lengthy 
algebra for t = r = const we obtain 



Ci-C: 



in 



T 



■,ni 



(73) 



Then performing numerical evaluations we arrive at in- 
equalities 

0< Cl 7 ClII <r^ (#) "^ (1) < 0.23 (74) 



0.15 



dT 

, Tdm 
mdT 



< T 



given gluoi 
dm (w) ~ JW) ( 1 ) ^ Ci - C: 



Cm 



for < T4tf < m , as in the given gluon model, and 



dT 



(1) 



< 



m 



C 



in 



(75) 



(67) 

for ^ < 0, as in case of the baryon-less matter described 
in the relativistic mean-field hadronic models (see [111, 
UH). Thus we see that in our gluon model values of 
C evaluated following expressions (|43]) and ([70]) deviate 
from each other by less than 23%. 



IV. RESULTS OF CALCULATIONS OF SHEAR 
AND BULK VISCOSITIES 

Analytical expressions for the shear viscosity coincide 
in all three ansatze provided one uses the same value of 
the averaged relaxation time f. 

For the bulk viscosity, the three expressions ([43]) , (|59l) , 
and (|70p coincide in the limit m — const but they differ in 
a general case. Also, results for £ obtained within ansatz 
I [Eqs. (03]) and g7j], ansatz II [Eqs. ([59]) and ([6T]) ]. 
and ansatz III [Eq. ([70p ] approximately coincide in the 
nonrelativistic limit. In order to estimate differences in 
the bulk viscosities calculated following Eqs. ([47]) . (|6ip . 
and (|70p . we further perform numerical evaluations. 

First, we should choose the value of the averaged re- 
laxation time f . In Ref. [l~6l | two parametrizations were 
used. The first parametrization, 



a g 2 Tln(2c/g 2 



(76) 
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with ao = Nr 1 47r ~ 0.2387 and a tuning parameter c 
(where in [16( we have chosen 2c = g 2 (T c )), is based 
on nonperturbative evaluations [s|. The deficiency of 
this parametrization is that it does not reproduce the 
appropriate perturbative limit for g 2 <C 1. The sec- 
ond parametrization, previously applied in Ref. flij . 
r _1 = a v /(32ir 2 ) g 4 T \a.{a ri 'K / g 2 ) with a tuning parame- 
ter a^, allows one to reproduce an appropriate perturba- 
tive limit for j 2 < 1. In recent papers [17| . the authors 
slightly modified the latter parametrization as 

f- 1 =ai 3 4 Tln( a2 /. 9 2 ), (77) 

introducing a tuning parameter <x 2 ; here ai = 0.02587 ~ 
a J) /(32 7 r 2 ). We use a 2 ==q 2 (7V) since for a 2 = (n* /T) 2 < 
g 2 {T c ) used in Refs. [15|, [l7| the result ([77)) becomes 
negative in some temperature interval above T c . 





V 
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■ 

• 
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FIG. 2: (Color online) The ratio of the gluon shear viscosity 
to the entropy density as a function of reduced temperature. 
Solid and dashed lines are computed with the averaged relax- 
ation time obeying Eq. (|77|) and Eq. (|76[) . respectively, with 
2c = a 2 = g 2 {T c ). The horizontal dotted line is the AdS/CFT 
bound rj/s = 1/4.7T [25|. The lattice gauge SU(3) data with 
16 3 x 8 and 24 3 x 8 lattice are from Ref. [f| (triangles and 
squares) and [f| (filled circles). The shaded band corresponds 
to the perturbative result of Ref. [2(|. The vertical dotted 
line shows the transition temperature. 

The temperature dependence of the rj/s and Q/s ratios 
is presented in Figs. [5] and [3J respectively. The lattice 
data were obtained using the lattice entropy from the 
old paper Q. We corrected these data for T/T c > 2 in 
accordance with the new lattice QCD result [H , which re- 
sulted in an increase of the lattice points in Figs. [5]and[3] 
by about 20% compared to the original paper @. As we 
have noted, all three ansatze result in the same values 
of the shear viscosity. The dashed and solid curves in 
Fig. [21 corresponding to the two parametrizations ([76) 
and (|77[) of the averaged relaxation time f , have simi- 
lar trends but different absolute values. With selected 
parameters c and a 2 in (|76j) and (|77[) the curves exhibit 
discontinuity at the critical temperature T c . Using (|77[) 



we can accurately describe the perturbative tail at very 
high temperatures. 

In Fig. [3j calculations using different expressions for 
the bulk viscosity are compared with the lattice result. 
One should note that although the presented lattice data 
are not very conclusive due to large error bars, never- 
theless, one may make some conclusions. The £/s ratios 




FIG. 3: (Color online) The ratio of the gluon bulk viscos- 
ity to the entropy density as a function of temperature. The 
solid and long-dashed curves are calculated following Eq. (|43[) 
(ansatz I) with relaxation time from (|77[) and (|76|l . respec- 
tively. The short-dashed curve is the ansatz II result [Eq. 
([59)) ] for the relaxation time from Eq. ([77)) . The dash-dotted 
line corresponds to the ansatz III calculations [Eq. (|70Jj with 
parametrization (|77[) . Lattice points are taken from @] (tri- 
angles) and 0] (filled circles). The shaded band corresponds 
to the perturbative result from Ref. [27fl . The vertical dotted 
line shows the transition temperature. 

calculated with the relaxation time (J77J within ansatze I, 
II, and III are presented by the solid, short-dashed, and 
dash-dotted curves, respectively. The results for ansatz 
I [Eq. d33)), solid curve] and ansatz III [Eq.fJTO]), dash- 
dotted curve] are very close to each other, in agreement 
with ([74"]) . whereas the ansatz II result [ ([5T))) . short-dashed 
curve] differs from them significantly and overestimates 
lattice data in the perturbative regime. The result ([5^)1 
of ansatz II with the relaxation time (J7BJ (see Fig. 4 of 
[lt| ) proves to be close to those presented by the solid and 
dash-dotted curves. On the other hand, the result (|4"3"[) 
(long-dashed curve, ansatz I) with the relaxation time 
(|76l) perhaps underestimates the lattice data, whereas in 
the nonperturbative regime (at T ~ T c ) the long-dashed 
curve using the relaxation time (|76p is closer to the data 
points than the solid curve using the relaxation time (|77p . 

V. CONCLUSIONS 

We obtained expressions for the gluon shear and bulk 
viscosities within the relaxation time approximation by 
making use of an averaged value for the relaxation time 
t = t = f. Under this assumption, the expressions for 
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the shear viscosity within ansatz I [Eq. (|35j) ]. ansatz 
II, and ansatz III [Eq. ([55] coincide, which indicates 
robustness of the shear viscosity to the ansatz reductions 
performed. 

In contrast, the expressions for the bulk viscosity, ([43]). 
(|59|) . and (f70|) . obtained within ansatze I, II, and III, re- 
spectively, significantly differ. Our numerical analysis 
demonstrates that the results (|4"3")l and ((70)) for gluons 
those mass increases with temperature, which use the 
Landau-Lifshitz condition (26), are close to each other 
in the whole range of temperatures whereas the result 
([55]) using a modified Landau-Lifshitz condition (27) de- 
viates from them significantly. One should stress that 
the result ([59]) is also recovered if in making variations 
one does not vary quantities which depend on the distri- 
bution function only implicitly, such as E within ansatz 
I. The later approximation is fully justified only for non- 
relativistic systems. 

Among the results for the bulk viscosity, Eq. (|4*3")l 
seems to us most theoretically established. Nevertheless, 
all evaluations of the shear and bulk viscosities obtained 
in the framework of the relaxation time approximation 
can be considered only as rough estimations. In order 
to perform more established calculations, one should go 
beyond the scope of the relaxation time approximation. 



However, such calculations are much more involved than 
estimations presented in the given work and have not yet 
been carried out for systems with strong interactions. 

Viscosity coefficients in a weakly coupled scalar field 
theory at arbitrary temperature can be evaluated directly 
from first principles without reference to the relaxation 
time approximation. This has been done by considering 
the expansion of the Kubo formulas in terms of ladder 
diagrams in the imaginary time formalism. In a theory 
with cubic and quartic interactions, the infinite class of 
diagrams which contribute to the leading weak-coupling 
behavior arc identified and summed. The resulting ex- 
pression is reduced to a single linear integral equation, 
which is shown to be identical to the corresponding re- 
sult obtained from a linearized Boltzmann equation sim- 
ilar to those which arise when the Boltzmann equation is 
treated in the relaxation time approximation, as was first 
noted in Refs. [28|,[29j]. Unfortunately, a similar analysis 
for more general cases is unavailable. 
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